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ABSTRACT 

We present a detailed study of the X-ray, optical and radio emission from the jet, lobes and core 
of the quasar PKS 2101-490 as revealed by new Chandra, HST and ATCA images. We extract the 
radio to X-ray spectral energy distributions from seven regions of the 13" jet, and model the jet X-ray 
emission in terms of Doppler beamed inverse Compton scattering of the cosmic microwave background 
(IC/CMB) for a jet in a state of equipartition between particle and magnetic field energy densities. 
This model implies that the jet remains highly relativistic hundreds of kpc from the nucleus, with 
a bulk Lorentz factor r ~ 6 and magnetic field of order 30 fiG. We detect an apparent radiative 
cooling break in the synchrotron spectrum of one of the jet knots, and are able to interpret this in 
terms of a standard one-zone continuous injection model, based on jet parameters derived from the 
IC/CMB model. However, we note apparent substructure in the bright optical knot in one of the 
HST bands. We confront the IC/CMB model with independent estimates of the jet power, and find 
that the IC/CMB model jet power is consistent with the independent estimates, provided that the 
minimum electron Lorentz factor 7 m i n > 50, and the knots are significantly longer than the jet width, 
as implied by de-projection of the observed knot lengths. 

Subject headings: galaxies: active - galaxies:jets - quasars individual (PKS 2101-490) 



1. INTRODUCTION 

The first Chandra observations of the quasar 
PKS 0637-752 revealed a bright X-ray jet extending 
12" from the quasar core (> 500 kpc de-projected, as- 
suming a jet viewing angle 9 < 9°, as evidence by 
the observed proper moti ons of the pc-sca le jet us- 
ing a modern cosmology (jLovell et al.1 I2006T )). associ- 
ated with the previously known radio jet, but with 
an unexpectedly high X-ray to radio flux density ratio 
(jSchwartz et al.1 120001: IChartas et al.l 120001) . Since then, 
tens of quasar jets have been found to possess X-ray jets 
with similarly high X-ray to radio flux density ratios (e.g . 
Harris fc Krawczvn ski 2002: iSambruna et alJl2002l 1 20041: 
Marshall et al.l 120051 1201 It iKataoka fc Stawarzl I2OO5I : 
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iMassaro etldl 120111) . The strong X -ray emission from 
kpc-scale quasar jets such as that of PKS 0637-752 is 
hard to explain in terms of standard emission mecha- 
nisms such as therm al bremsstrahlung or synchrotron self 
Comp ton emission (|Chartas et al.l 120001 : iSchwartz et~aTl 
2000) . A popular explanation for the strong X-ray emis- 
sion is the beamed, equipar tition IC/CMB model pro- 
posed for PKS 0637 -752 by iTavecchio etaLl (|2000l ) and 
iCelotti et al.l (|2001l ). in which the flow velocity is as- 
sumed to be highly relativistic and directed close to the 
line of sight. A relativistic jet velocity increases the en- 
ergy density of the cosmic microwave background (CMB) 
in the rest frame of the jet plasma, thereby increasing the 
X-ray emissivity produced via inverse Compton scatter- 
ing of CMB photons. The small jet viewing angle im- 
plies that the emission is Doppler beamed towards the 
observer. The appeal of this model is largely due to 
its simplicity and consistency with equipartition between 
magnetic and particle energy densities in the emitting 
plasma. From here on, we refer to the beamed, equipar- 
tition IC/CMB model as simply the IC/CMB model. 

A number of uncertainties and potential problems for 
the IC/CMB model have been identified: (1) There is 
no conclusive theoretical or empirical justification for 
the assumption of equipartition in jet plasma, although, 
for a given jet speed, the equipartition condition min- 
imises the plasma energy density and jet power. (2) 
The IC/CMB model requires jet Lorentz factors of or- 
der r ~ 5 — 25 on scales of hundreds of kpc from 
the core. Such large jet Lorentz factors are inconsis- 
tent with studies of the radio e mission from large-scale 
jets and counter-jets (r 5 , 1.5 iWardle fc Aaronl 119971; 
iMullin fc Hardcastlel l2009f ) ■ A suggested solution to 
this problem invokes velocity structure across the jet — 
the so-called "spine-sheath" model, in which the radio 
emission from jets in FRII radio galaxies is dominated 
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by a slower moving sheath, whilst the emission from 
quasar jets with small viewing angles is said to be domi- 
nated by the Dop pler boosted radiation from a fast mov- 
ing sp ine (see e.g. lHardc astle 200(| iMullin fc: Hardcastld 
2009). (3) It has been argued that, due to the long 
cooling timescale for the 7 ~ 100 electrons responsible 
for the IC/CMB X-ray emission, radiative and adiabatic 
losses alone cannot account for the r apid drop in X- 
ray s u rface brightnes s outside the knots (ITavecchio et all 
120031: IStawarz et all 12004 ISiemiginowska et al.l I2007D 
(4) The IC/CMB model makes strong, testable predic- 
tions about the redshift dependence of kpc-scale X-ray 
jets. Specifically, the model predicts that the X-ray sur- 
face brightness should be redshift independent, because 
the CMB energy density increases as (1+z) 4 which bal- 
ances the usual (1+z) -4 decrease of surface brightness 
with redshift. Therefore, the ratio of X-ray to radio 
surfac e brightness should be a strong fun ction of red- 
shift (lSchwartzll20ol Marshall et al.l l20lTh . So far, the 
search for th e predicted red s hift dependence has been 
unsu ccessful ([Marshall et al.1 120111 : iKataoka fc Stawarz) 
20051) . (5) It has been argued (e.g. lAtovan fc Dermerl 
2004) that the very large jet powers d erived from the 
IC/CM B model (> 10 48 ergs/s, e.g. ITavecchio etafl 
(2000)) are prohibitively large. Such high jet power is 
disfavoured because 10 48 ergs/s is equal to the Edding- 
ton luminosity of a 10 10 M Q black hole, and such high jet 
power is an outlier when c ompared to samples of FRII 
radio galaxies such as the iRawlings fc Saunders! ()1991[ ) 
sample, for which the largest estimated jet power is of 
order 10 47 ergs/s. (6) A number of sources show a de- 
creasing X-ray to radio flux density ratio along the jet, 
which, if the IC/CMB model is correct, implies decel- 
eration must be taking place on h undreds of kpc-scales 
(jGeorganopoulos &: Kazanasl 12004) . However, it is not 
clear how the gradual deceleration can occur, and there 
is no evidence for deceleration in lobe-dominated ra- 
dio galaxies, whi ch may be expected in such a model 
(|Hardcastlell2006l) . 

Despite the numerous concerns surrounding the 
IC/CMB model, none of the issues listed above is cur- 
rently seen as definitively refuting the model, and it con- 
tinues to receive attention in the literature as the likely 
candidate for the quasar jet X-ray emission mechanism. 
In this paper, we critically assess the application of the 
IC/CMB model to jet X-ray emission of PKS 2101-490. 

PKS 2101-490 was fir st repo r ted as a bright flat spec- 
trum radio source by lEkersI (|1969D . IMarshall et"aT1 
(2005) reported a redshift of z sa 1.04 for this source, 
based on an unpublished Magellan spectrum (see also 
the discussion in <J6]). The spectroscopic redshift deter- 
mination is robust, with uncertainty of approximately 
±0.003. Further details of the spectroscopic observations 
and data analysis will be presented in an upcoming paper 
(Gelbord & Marshall, in prep.). 

Studies at the Australia Telescope Compact Array 
(ATCA) re vealed signifi cant radio emission on arcsec- 
ond scales (|LoveU 119971) . For this reason, PKS 2101- 
490 was included in the Chandra snapshot survey of 
flat spec trum radio quasars wit h arcsecond sc a le ra - 
dio jets (|Marshall et all 120051 ). IMarshall etldl (|2005l ) 
presented an 8.6 GHz ATCA image along with a 
5 kilosecond snapshot Chandra image (Chandra Ob- 



TABLE 1 
Observation Information 



Instrument 


Band 


Mode 


Date 


ATCA 


4.8 GHz 


1.5A/6A 


May 25/Sep 2 2000 


ATCA 


8.64 GHz 


1.5A/6A 


May 25/Sep 2 2000 


ATCA 


17.73 GHz 


6C 


May 10 2004 


ATCA 


20.16 GHz 


6C 


May 10 2004 


HST 


F814W 


ACS/WFC 


Mar 8 2005 


HST 


F475W 


ACS/WFC 


Mar 8 2005 


Chandra 


0.5 - 7 keV 


ACIS-S3 


Dec 17 2004 



sID 3126) that revealed significant X-ray emission as- 
sociated with the 13" eastern radio jet. Based on 
the results of the snapshot survey, and its morpho- 
logical similarity to PKS 0637-752, PKS 2101-490 
was selected along with a number of other sources : 
PKS 1421-490 ([Godfrey et al|[2TMlGeTbord et all2005j) : 
PKS 1055+201 (ISchwartz et al.l l2006bD: PKS 0208^512 
(iPerlman et all [2011f); P KS 1202-262 (IPerlman et all 
I2011J); PKS 0920-397 ([Schwartz et al.l 120101) : and 
PKS 1030-357, as a target for deeper follow-up obser- 
vations with Chandra, HST and ATCA. Here we present 
a detailed physical analysis of the jet, hotspot and lobes 
of this source based on new ATCA, Chandra and HST 
follow-up images. 

In S}2]we describe the observations and data reduction. 
In we discuss the method and results of modelling 
the spectral energy distributions of spatially separated 
jet knots in terms of the IC/CMB model. In SQ]we com- 
pare independent jet power estimates with that obtained 
from the IC/CMB model of jet X-ray emission. In §5\ 
we discuss the optical emission detected from one of the 
jet knots, and present an interpretation of this in terms 
of a broken power law synchrotron spectrum. In SjS] we 
discuss the X-ray spectrum of the quasar core. In Sj7]we 
discuss the radio and X-ray emission from the lobes, in 
particular, we model the lobe SEDs in terms of inverse 
Compton scattering of the CMB. In $8] we present the 
conclusions and final remarks. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Overview 

Figure Q] illustrates the radio structure of the source 
and the naming convention used for the various features 
in the radio maps. Figure [2] illustrates the polarization 
characteristics of the jet. We have extracted radio, opti- 
cal and X-ray flux densities from the seven major emis- 
sion regions identified in Figure [TJ In this section we 
describe the observations in each waveband, as well as 
the methods used to measure flux densities and sizes for 
the individual jet knots. Table Q] lists the observational 
information for all data used in this study. 

2.2. Radio 

PKS 2101-490 was observed with the ATCA at 
4.8 GHz and 8.64 GHz in two configurations, 1.5A and 
6A, on May 25 2000 and September 2 2000 respec- 
tively, and in a single configuration (6C) at 17.7 GHz 
and 20.2 GHz on May 10 2004. In each case, a full 12 
hour synthesis was obtained, recording 128MHz band- 
width in all four polarization products. Regular scans on 
the nearby phase calibrator PKS 2106-413 were sched- 
uled throughout the observations, as well as scans on 
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Fig. 1.— ATCA images at 4.8 GHz (left) and 17.7 GHz (right) showing the radio structure of the jet and lobes of PKS 2101-490. Also 
shown is the naming convention used for various components of the jet. The scale of this image is 8.1 kpc/". Contours are separated by a 
factor of 2 in surface brightness. In the 4.8 GHz image, the lowest contour is 0.42 mjy/beam and the beam FWHM is 2''24 X 1"84. In the 
17.7 GHz image, the lowest contour is 0.15 mjy/beam and the beam FWHM is 0"79 X 0'.'54. 




Relative Right Ascension (arcsec) 

Fig. 2. — Same as Figure [l] (right) with polarization E-vectors overlaid. The E-vectors are scaled such that a length of 1 arcsecond 
corresponds to 1 mjy/beam polarized surface brightness. 



the ATCA flux calibrator PKS 1934-638. Standard cal- 
ibration and editing procedures were carried out using 
the MIRIAD data analysis package. Following calibra- 
tion, the data were exported to DIFMAP and several 
imaging/self-calibration iterations were performed. The 
data were both phase and amplitude self-calibrated. 

2.2.1. Radio Knot Flux Density and Size Meausrements 

The spectrum of the entire jet (excluding lobe emis- 
sion) between 4.8 GHz and 20.2 GHz is well described 
by a power law with spectral index a = 0.81 ± 0.01 (flux 
density S u oc v~ a ). Each of the four flux density mea- 
surements are within 1% of the best fit power law (see 
Table [2]), giving us confidence in the flux density and 
spectral index measurements for the entire jet. However, 
inspection of individual knot spectra indicates that the 
systematic uncertainties in flux density measurements for 
individual knots are significantly greater than the off- 
source RMS. There are a number of factors contribut- 
ing to the systematic uncertainty in the measurement of 
the relative strengths of individual knots, including the 
uniqueness problem that results from gaps in the (u, v)- 
coverage due to the small number of telescopes (see e.g. 
iWalkerl [1995D . and weak inter- knot emission that is de- 
tected at different levels in each band, as a result of the 



different (u, v)-coverage at each frequency. Also, some of 
the knots at 8.6 GHz are only marginally resolved, so that 
flux density measurements of individual knots becomes 
difficult. These issues, and estimat ion of systemat ic un- 
certainty are discussed in detail in IGodfrevl ([2009D . We 
are therefore unable to determine accurately the spectral 
index for each individual knot, and are instead forced to 
assume that each knot has the same spectral index as 
that of the entire jet, which has been determined accu- 
rately. 

Knot sizes for jet components Kl - K4 were measured 
by fitting 2D Gaussian models to jet knots in the 17.7 
and 20.2 GHz maps, using the AIPS task JMFIT. Knot 
6 and the hotspot are partially blended even at 20 GHz, 
and for this reason we determined their flux density and 
size at 17.7 and 20.2 GHz by fitting elliptical Gaussian 
models to these components in the (u,v)-plane, using the 
modelfit function of DIFMAP. The knot size estimates 
for all knots are based only on the measurements from 
the highest resolution maps (17.7 and 20.2 GHz). At 8.6 
GHz, for knots 2, 3 and 4, flux densities were obtained by 
integrating the surface brightness within a region encom- 
passing the knot emission. For knot 6 and the hotspot, 
the flux densities at 8.6 GHz were determined by fitting 
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TABLE 2 

Characteristics of Spatially Resolved Jet Knots and Lobes 



Flux Densities De-convolved Dimensions 



Knot ID F 4 . 8 GHz Fs.6 GHz F17.7 GHz F 2 Q.2 GHz Fs06nm F475 nm F x kcV </>Maj a 0Min a Vol b 

mjy nijy mjy mjy njy njy njy mas mas cm 3 



Core 


720 ± 40 


830 ± 40 


500 ± 50 


470 ± 50 


2.32 x 10 5 


1.59 x 10 5 


87 ±3 










Knot 1 






5.1 ± 0.5 


4.7 ±0.5 


60 ± 10 


30 ± 10 


< 0.2 


400 ± 100 


250 ± 10 


3 x 


10 65 


Knot 2 




4.6 ±0.5 


2.3 ±0.2 


2.2 ± 0.2 


< 30 


< 30 


< 0.15 


400 ± 100 


<150 


1 X 


10 65 


Knot 3 




5.1 ±0.5 


3.5 ± 0.3 


2.8 ±0.3 


< 30 


< 30 


0.5 ±0.15 


550 ± 100 


330 ± 100 


7 x 


10 65 


Knot 4 




9± 1 


4.7 ±0.5 


4.2 ± 0.4 


< 50 


< 50 


1.3±0.3 


2500 ± 500 


300 ± 50 


3 x 


10 66 


Knot 5 






1.0 ±0.2 


0.8 ±0.2 


< 30 


< 30 


0.2 ±0.1 


400 ± 100 


300 ± 100 


4 x 


10 65 


Knot 6 




11 ± 1 


6.5 ± 0.6 


6.3 ±0.6 


180 ± 10 


90 ± 10 


0.75 ±0.2 


600 ± 50 


400 ± 20 


1 X 


10 66 


hotspot 




30 ±3 


16 ± 1.5 


14 ± 1.5 


< 30 


< 30 


< 0.16 


340 ± 30 


240 ± 20 


2 x 


10 65 


Entire Jet 


115 ±1 


71 ± 1 


40 ± 0.5 


36 ±0.5 






3.3 ±0.4 










Lobe 


16 ±2 


8± 1 


4.5 ± 0.4 


4 ±0.4 






0.5 ±0.2 


8600 


6500 


3 x 


10 69 


Counter-Lobe 


65 ±6 


33 ±3 


17 ±2 


14 ± 1 






1.5 ±0.2 


12000 


9400 


9 x 


10 69 



These arc the FWHM Gauasiai 
was performed in the (u,v)-planc 
quoted sizes arc simply the size c 
k These arc apparent (projected 
diameter of the ellips' 
to the line of sight. 



1 knot sizes, based on Gaussian fits to the ra< 
(sec i|2.2.1t . For all other jet knots, a Gaussi 
f the flux extraction regions — no fitting was 
I volumes calculated assuming ellipsoidal gco: 
If the knots arc associated with stationary features ii 



io data (except in the 
in model was fit to th 
done for the lobes, 
ictry, i.e. V = (tt/6) 
the jet, the de-projected 



of the lobes). In the 



eq 



ere D e q is the cquatori 
'ill be greater by a facto 



nd the hotspot, the Gauss 



.1 diame- 
■ of (1/f 



elliptical Gaussian models to these components in the 
(u,v)-plane. Due to the limited resolution at 8.6 Ghz, we 
were unable to accurately determine the flux density for 
knots 1 and 5. 

The size estimates from the 17.7 GHz and 20.2 GHz 
images were consistent for some knots. However, in some 
cases the knot dimensions differed by as much as a factor 
of 2. This is another indication that image fidelity is 
questionable in the high frequency maps, and systematic 
errors resulting from the deconvolution/self-calibration 
process are significant. The knot sizes listed in Table 
[2] are averages of the parameters determined from the 
17.7 GHz and 20.2 GHz images. The uncertainties in 
knot size are taken to be half the difference between the 
17.7 and 20.2 GHz measurements. 

2.3. Hubble Space Telescope Optical Observations 

PKS 2101-490 was observed with the Advanced Cam- 
era for Surveys (ACS) Wide Field Camera 1 (WFC1) 
on the Hubble Space Telescope in two filters (F475W 
and F814W) on March 8 2005. A total exposure time of 
2.3 kiloscconds was obtained in each filter. A sub-pixel 
dithering pattern, with CR-SPLIT images at each of 3 
positions along a chip diagonal, was utilized to eliminate 
bad pixels and allow us to maximize the angular resolu- 
tion, as the ACS/WFC does not fully sample the PSF 
at either 4750 or 8000 A. An ORIENT was chosen such 
that the jet did not fall within 25 degrees of a diffrac- 
tion spike. The data were retrieved from the Multi- 
mission Archive at Space Telescope (MAST) websitcFl. 
however, multiple peaks in the pipeline drizzled image 
of the quasar core suggested that a re-reduction of the 
data was required. In addition, the pipeline does not take 
into account sub-pixel dithering, which was performed to 
in order to recover PSF information from the undersam- 
pling of the HST PSF. The re-reduction involved running 
MULTIDRIZZLE with a smaller PIXFRAC and scale, 

13 HST archive http://archive.stsci.edu/ 



so that the images could be sub-sampled to 0.0247 arc- 
second pixels, and checking the alignment of the images 
with TWEAKSHIFTS. The best reference files indicated 
in the HST archive were used for the re-reduction. The 
position of the optical quasar core was aligned with the 
position of the radio core. This required shifting the op- 
tical data by approximately 0'.'4. 

We used TinyTIM simulations (IKrist fc Burrowsl[l994l: 
ISuchkov fc Kristl 119981: IKrist fc Hookl I2004D toobtain 
subsampled PSF simulations for both bands. For those 
simulations we assumed an optical spectrum of the form 
F v oc i/ -1 , although experience has shown that the PSF 
shape is not heavily dependent on spectral slope. Ro- 
tating the PSF to a north-up position for use with the 
drizzled images required independently rotating the two 
axes as the WFC detector's axes are not completely or- 
thogonal on the sky. We normalized the PSF to the total 
number of counts in a 2" circular aperture centred on 
the source. This takes advantage of the fact that charge 
"bleed" on the AC S is linear and c harge is conserved for 
a saturated source (Gilliland 2004), and enables optimal 
matching of the PSF's outer portion to what is observed. 
Because the quasar core was saturated in the individual 
exposures, this inevitably led to negatives in the central 
pixels when PSF subtraction was done; however, given 
the small residuals in other places plus the smooth off- jet 
isophotes, we believe the result is reliable. Subtraction 
of the PSF allowed us to look for optical jet components 
within 1-2", and resulted in a detection of Knot 1 in both 
bands. 

Figure [3] shows the resulting optical maps with radio 
contours overlaid. In these images, the HST data have 
been smoothed with a 0.3 arcsecond FWHM Gaussian to 
better show the optical emission from knot K6, which is 
clearly detected in both filters. It is interesting to note 
that the F814W image of K6 reveals an elongated, double 
structure that is not apparent in the F475W image of K6. 
We further note that the position of the optical emission 
from K6 is coincident with the radio position to within 
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TABLE 3 

Knot 6 (K6) de-convolved sizes 



Band 


Cross-jet a 


Along-jet a 




(arcsec) 


(arcsec) 


ATCA 20.2 GHz 


0.4 ±0.02 


0.6 ±0.05 


HST F814W 


24+ 03 
u ' z4 -0.24 


0.4 ±0.1 


HST F475W 


24+ ' 03 
u ' z4 -0.24 


n O+0.03 
u - z -0.2 


Chandra 0.5-7 kcV 


< 0.6 


< 0.6 



Fig. 3.— Hubble Space Telescope images of PKS 2101-490 (grayscales: Top = F814W; Bottom = F475W) with 20.2 GHz ATCA radio 
contours overlaid. The HST images have been sub-sampled with pixel size 070247, PSF subtracted from the quasar core, and smoothed 
with a 0''3 FWHM Gaussian to better show the optical counterpart to knot K6. The diffraction spikes are merely the residuals after PSF 
subtraction from the quasar core (see text). 

the uncertainties associated with the optical-radio image 
alignment. We consider the interpretation of the optical 
data in !}5l 



2.3.1. Optical Flux Density Measurements 

Knots 1 and 6 are the only jet features detected 
in the optical images. Optical flux densities were 
measured using standard aperture photometry tech- 
nique s. The approp ri ate ap erture corrections were taken 
from ISirianni et all (|2005T ) Table 3, and the appropri- 
ate ex tinction corrections were taken from ISirianni et al.l 
(|2005h Table 14 assuming E(B-V) = 0.039 at the position 
of PKS 2101-4903. The optical flux densities and upper 
limits are given in Table [2j 

2.3.2. The extent of optical emission from Knot 6 

We estimated the size of the optical emission from K6 
as follows: first, the optical images were convolved with 
0715 FWHM Gaussian. We then produced integrated 
profiles both parallel and perpendicular to the jet axis, 
for the knot K6, a nearby star, and the Tiny Tim gener- 
ated PSF using the ds9 projection capability. From the 



These are the intrinsic FWHM Gaussian knot 

ZCS, Obtained Via 6>K6, intrinsic — y®l 



Br- 



J K6 — ^PSF' 

except in the case of F814W along-jet, which ex- 
hibits a double peak structure. In that case, the 
extent of the emission region was estimated by in- 
spection of the longitudinal knot profile (Figure^. 

integrated profiles, we measured the FWHM of eactp 
feature to be Ok6, ©star and Opsf- The profiles of the 
star and TinyTim PSF are in good agreement with each 
other, and indicate that in the smoothed maps, the PSF 
FWHM is Bpsf ~ Sstar ~ 07 18 for both the F475W 
and F814W bands. We then calculated the de-convolved 
(intrinsic) extent of the optical emission associated with 

the knot K6 as 6 K6, intrinsic = 



©K6 



e 2 

W PSF/Sta 



The 



This value for E(B-V) was obtained using the 
NASA Extragalactic Database extinct ion ca lculator 
http: / /nedwww. ipac . caltech.edu/forms/calculator . html 
which is based on the Galactic reddening maps of Schle gel et al.l 



measured PSF sizes in the smoothed HST maps corre- 

15 The longitudinal profile in the F814W band appears double 
peaked, so in that case, we estimated the extent of the optical 
emission region simply via inspection of the profile (Figure [4}. 
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spond to an un-smoothed PSF of O'/ll in both bands, as 
expected for ACS/WFC. We note that the profile of K6 
along the jet in the F814W band appears double peaked, 
in contrast to the F475W profile which is single peaked 
(see Figures [3] and 0]) . 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 
Relative Distance Along Jet [arcsec] 



Fig. 4. — Longitudinal jet profiles of knot 6 in both HST bands 
along a position angle of 150°, integrated across a region of width 
0''5. These plots serve to illustrate the difference in the apparent 
knot length in each band, and the relative alignment of the peaks. 
Note that the peak in the F475W profile lies between the peaks in 
the F814W profile. The HST maps were smoothed by a Gaussian 
with CK'15 FWHM before producing these profiles. The profile of 
the PSF, smoothed by a Gaussian with 0'.'15 FWHM, is represented 
by the dashed curve in each plot. The error bars are based on the 
rms of profiles made in off-source regions. 

The optical knot appears smaller than the associated 
radio knot: in both bands the cross-jet width of the opti- 
cal emission associated with K6 is significantly less than 
the measured width in the radio band (Table [3J . The 
length of K6 parallel to the jet in the F814W map ap- 
pears marginally resolved, while the length of K6 in the 
F814W map is clearly resolved, and appears to consist of 
two peaks (Figured]). A discussion of the optical emis- 
sion from K6 is presented in Sj5] 

2.4. Chandra X-ray Observations 

PKS 2101-490 was observed with the Chandra X-ray 
observatory on 17 December 2004 (Cycle 6) using the Ad- 
vanced CCD Imaging Spectrometer (ACIS-S) for a total 
exposure time of 42 kiloseconds (Chandra ObsID 5731). 
To reduce the effect of pile- up in the quasar core, a 1/4 
subarray mode was used with a single CCD, so that the 
frame time was 0.8 seconds. The source was positioned 
close to the readout edge of the CCD to reduce the ef- 
fect of charge transfer inefficiency A new Event 2 data 



file was made with pixel randomization removed, and the 
X-ray and radio core positions were aligned, requiring a 
shift of the X-ray data by approximately 0'.'4. The data 
were restricted to the energy range 0.5 - 7 keV. 

2.4.1. X-ray Flux Density Measurements 

The X-ray flux densities for individual knots were ob- 
tained by calculating the background-subtracted counts 
in each knot region and multiplying by the conversion 
factor. The conversion factor was estimated by fitting 
an absorbed power law spectrum to the X-ray counts ex- 
tracted from the whole jet. The X-ray spectra of the jet 
and lobe were fit using the Sherpa software package by 
minimising the Cash statistic (related to the log of the 
likelihood). The instrument response functions (RMF 
and ARF) were determined from the CALDB calibra- 
tion database appropriate for the position of the source 
on the ACIS-S3 chip. In all model fits, the neutral Hy- 
drogen column density was fixed at the Galactic value 

3.4 x 10 20 cm -2 as determined from the COLDER col- 
umn density calculator provided by the Chandra X-ray 
Center. The uncertainties for the flux density and spec- 
tral index were calculated using the "Covariance" rou- 
tine in Sherpa. The results of this routine are valid pro- 
vided that the surface of Log Likelihood is approximately 
shaped like a paraboloid. The "Interval-Projection" rou- 
tine in Sherpa was used to verify that this condition was 
met. 

We extracted a total of 138 counts from the jet having 
energies in the range 0.5 keV to 7.0 keV. Fitting a power 
law gives 1 keV flux density S = 3.3±0.4 nJy and spectral 
index ax = 0.85±0.2. This implies the conversion factor 
is G « 1.0 /jJy/count/s for a spectral index of a « 0.85. 
In the counter lobe we extracted a total of 55 counts 
in the energy range 0.5 keV to 7.0 keV. Following the 
same procedure as for the jet, we find flux density S — 

1.5 ± 0.2 nJy at 1 keV and spectral index a — 1.3 ± 0.3, 
implying the conversion factor is G = 1.15 /iJy/count/s 
for a spectral index a ~ 1.3. These estimates of the 
conversion factor are consistent with the predictions of 
the Chandra proposal planning toolkit. 

Figure [SJ^. is a comparison of X-ray and radio structure 
in the jet of PKS 2101-490 and Figure[5)3 illustrates the 
regions used to extract X-ray counts for the knots. The 
regions shown in Figure 03 (except for the hotspot ex- 
traction region — this region is discussed further below) 
have sides > l'/6. Knots 5, 6 and the hotspot are difficult 
to separate in the Chandra image. In order to estimate 
the counts associated with the hotspot and avoid con- 
tamination from knot 6, a region encompassing only one 
side of the hotspot is used (the side furthest from knot 
6). Due to the background in the vicinity of the hotspot 
and the possibility of contamination from knot 6, the few 
counts within this aperture may not be associated with 
the hotspot, and therefore the flux density estimate for 
the hotspot is an upper limit. An aperture correction 
of 2 is used when calculating the upper limit, since the 
flux extraction region encompasses only one side of the 
hotspot. 

In the Chandra image, knot 1 is blended with the wings 

16 The COLDEN column density calculator, available 
at http: //exc . harvard. edu/toolkit/ colden. jsp based on 
IDickev L Lockmanl lfl990i ). 
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Fig. 5. — (A) X-ray image binned to 0"0492 pixels and smoothed with a Of.' 7 FWHM Gaussian with 17.7 GHz ATCA contours overlaid. 
Colour scale is logarithmic between 0.0045 and 0.15 counts/beam. (B) Raw X-ray counts image binned to half a Chandra pixel width 
(0'.'246) with 17.7 GHz radio contours overlaid. Colour scale is logarithmic between and 25 counts. The black outlines mark the flux 
extraction regions used for X-ray flux density calculations. (C) X-ray image binned to 0'.'246 pixels and smoothed with Of.' 7 FWHM 
Gaussian. Colour scale is logarithmic between 0.01 and 0.18 counts/beam. Blue contours are from the 4.8 GHz ATCA image. (D) Raw 
X-ray counts image with 4.8 GHz radio contours overlaid. Colour scale is logarithmic between and 30 counts. The black outlines mark 
the flux extraction regions used to calculate X-ray flux densities associated with the lobes. 



of the X-ray core. To place a limit on the X-ray flux 
density from knot 1 , a sector of an annulus centred on the 
core was used, as shown in Figure [5j3. The background 
was estimated using the section of the annulus excluding 
knot 1. The data are consistent with zero counts from 
knot 1. 

2.4.2. The extent of X-ray emission from Knot 6 

The X-ray knot siz e is of great significanc e to models 
of jet X-ray emission (|Tavecchio et al.ll2003f) . Only knot 
6 has sufficient counts and is sufficiently isolated from 
other strong knots to allow an accurate estimate of the 
knot size. Using the CIAO task dmlist, events were ex- 
tracted from a 0'.'8 radius circular aperture centered on 
knot 6. This region was chosen to be large enough to 
include a large fraction of the counts from knot 6, but 
small enough to avoid contamination from neighbouring 
regions of the jet. Note that the encircled energy fraction 
within a circular aperture of radius 0'.'8 is > 85% for a 
spectral index of a ~ 0.8. A total of 29 events were ex- 
tracted from this region. We assume, for simplicity, that 
the knot surface brightness profile and the inner 0'.'8 of 
the Chandra PSF are both approximately Gaussian with 
variance cr^ 6 and crp gF Q // g respectively. We estimated 

the variance (a 2 ,, ) of the Chandra PSF within an 

v PSF,0. 8 y 

O'.'S radius aperture by extracting events within an 0'.'8 



circular aperture centered on the core, and calculating 
the standard deviation of event coordinates. Doing so, 
we find er pSF " 8 = 0'.'29 ± C/10. This comparison be- 
tween core and jet PSF is valid since the core and jet 
X-ray spectral indices are similar, and the core is not 
significantly affected by pile-up. We find that knot 6 is 
unresolved in both the jet and cross-jet directions. 

To obtain an upper limit to the size of knot 6 we use 
the st andard deviation distribution (|Kennedv fc Keeping) 
119511 ). We find that the 99% upper limit to the 



X-ray knot size is er b s 



0'.'4. The upper 



limit to the de-convolved (intrinsic) size of knot 6 
(taken as the Gaussian FWHM) is then Dkq < 

, / s 1/2 

2V21n2 a 2 h „ nm „ - a 2 „ = Of! 6 in both the jet 

^ obs, uppci PSF, 0. 8/ 

and cross-jet directions. Thus, the extent of the X-ray 
knot emission is less than or equal to the extent of the 
radio emission of knot 6. 

3. RESULTS 

In this section we present a comparison of the radio and 
X-ray images; we describe the methods used to model the 
radio to X-ray SEDs of individual knots; and we present 
the results of spectral modeling. 

3.1. Radio/X-ray longitudinal jet profiles 
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arcseconds from core 

Fig. 6. — Top: X-ray image binned with 1/10 Chandra pixel width (CK'0492) and smoothed with a Of.' 7 FWHM Gaussian to emphasize 
the X-ray jet structure. 17.7 GHz ATCA contours are overlaid. Also shown are the projection regions used to obtain the radio and X-ray 
jet profiles. Bottom: X-ray (dashed line) and radio (solid line) longitudinal jet profiles. These curves show the jet brightness integrated 
across the jet as a function of distance from the core along the jet ridge-line. Note that in this plot, the X-ray resolution (l'.'l) and radio 
resolution (0'.'65) are not equal. These profiles simply serve to illustrate the differences in jet structure at each wavelength. 



Figure [5] illustrates the radio and X-ray longitudinal jet 
profiles. The X-ray surface brightness is greatest near the 
end of the jet. This is in contrast to the trend that is 
seen m sources s uch a s 3C273 (e.g. iMarshall et alll200ll: 
Sambruna et all 120011). 0827+243 (IJorstad & Marscherl 
2004]) and 1127-145 (jsiemiginowska et alj 1200% 120071 ). 
where the X-ray brightness peaks closer to the core and 
the brightne ss of the knots decreases with dista nce from 
the core (see iGeorganopoulos fc Kaza nas 2004, for a dis- 
cussion of this phenomenon). 

3.2. Modeling the Spectral Energy Distributions of 
Spatially Separated Knots 

The spectral energy distributions of the jet knots 
are typical of quasar X-ray jet knots such as those in 
PKS 0637-752: A single or broken power law is un- 
able to fit the entire radio to X-ray SED, so that two 
spectral components are required to model the data (see 
Figure [7]). The X-ray spectral index of the entire jet 
a 05kcV = 0.85 ±0.2 is consistent with the radio spectral 
index of the entire jet an = 0.81 ± 0.01. The data are 
therefore consistent with an inverse Compton interpreta- 
tion for the X-ray emission. 

3.2.1. Synchrotron Self Compton Model 

As with other quasar X-ray jets (e.g. Schwartz et al.1 
2000), synchrotron self Compton (SSC) models for 
the X-ray jet emission in PKS2101-490 require sub- 
equipartition magnetic field strengths. The magnetic 
field strengths in the knots derived from SSC modeling 
are a factor of ^50 below the equipartition (minimum 
energy) values (|Godfrevl I2009D . Such sub-equipartition 



magnetic fields are generally disfavoured, because in that 
case, the internal energy, internal pressure and jet power 
(for a fixed jet speed) are orders of magnitude greater 
than in the case of equipartition magnetic fields. We 
note that if the jet is Doppler beamed, then the ratio 
of the equipartition magnetic field strength to the SSC 
derived magnetic field strength increases approximately 
in proportion to <5 2 /( 1+Q ). 

3.2.2. IC/CMB Modeling 

In this section we model the knot X-ray emission in 
terms of inverse Compton scattering of CMB photons in 
a highly relativistic jet directed close to t he line of sight 
(|Tavecchio et al.ll2000l: FCelotti et al.ll200l . 

Assumed form of the electron energy distribution: — In 
order to apply the analytic one- zone I C/CMB model 
(IDermerl 119951: Harris k, Krawczvnskil I2002t IWorralll 
120091 ) . we assume a single power law electron energy dis- 
tribution iV(7) = K e ^/~ a between a minimum and max- 
imum Lorentz factor, 7 m j n and 7 ma x, with a — 2.6 (cor- 
responding to synchrotron spectral index a = 0.8). We 
assume 7 max = 10 5 for all knots. The assumed value 
for 7 max allows the synchrotron spectrum to cut-off at 
a frequency below that of the HST observing frequen- 
cies. The value of the high energy cut-off is not well 
constrained, however, the conclusions drawn from this 
model are not sensitiv e to the assumed value of 7 max (see 
ISchwartz et al.|[2006al Appendix A). In contrast to the 
single power law electron energy distribution assumed 
here, in $5] we model the synchrotron spectrum of Knot 
6 in terms of a broken power law. We note that the re- 
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TABLE 4 
IC/CMB Modeling Results 



Knot ID 


F lkcV 
F 17.7 GHz 

[xlO" 8 ] 


<5 


Ben 
[ M G] 


[xlO -6 cm" 3 ] 


r e 

[xlO 46 erg 


;s/s] 


[xlO 46 ergs/s] 


Knot 1 


< 4 


< 6 


> 40 


> 0.6 


< 0.5 




< 3 


Knot 2 


< 6.5 


< 7 


> 40 


> 0.5 


< 0.2 




< 1 


Knot 3 


14 ±4 


7 


30 


0.2 


0.5 




3 


Knot 4 


28 ± 7 


7 


20 


0.1 


0.2 




1 


Knot 5 


20 ± 10 


7 


20 


0.2 


0.3 




1 


Knot 6 


11 ±3 


6 


30 


0.3 


0.7 




i 


Note. — 


Assumed model parameters: 


7min — 50, 7 ma3 


, = 10 5 , a = 


0.8, 





Jet viewing angle — 9 degrees (for dc-projection of knot volumes given in Table [2t , 
ratio of proton to electron energy densities — << 1. 
To calculate jet powers we have assumed T — 8. 



suits of the IC/CMB model are insensitive to the details 
of the high energy end of the electron energy distribu- 
tion, provided the break in the distribution occurs at 
7brcak >> 7min- Therefore, the assumption of a single 
power law distribution in this section does not affect the 
conclusions drawn from this model for Knot 6, nor does 
it alter the conclusions drawn in later sections. 

The low energy cutoff, 7minj is constrained so that 
the model does not produce optical IC/CMB emission 
above the HST upper limits or detections. This con- 
straint on 7 m i n is possible because an extrapolation of 
the IC/CMB spectrum from X-ray to optical frequencies 
lies above the optical upper limits, and in the case of 
Knot 6, is comparable to the detection level but with an 
incorrect spectral index. In addition to this constraint, 
7min must not be so high that there is no significant X- 
ray emission produced at 0.5 keV. As we show in the fol- 
lowing sections, the IC/CMB model requires a Doppler 
factor S ~ 6 at least in some parts of the jet. Assuming 

"m(n mh ~ L6 x 1011 x £ 2 7mm Hz ( see e.g. IWorral]|[2009h . 
we are able to constrain the value for 7m to lie in the ap- 
roxim ate range 10 < 7, n i n < 200. Muel ler fc Schwarta 



2009) find marginal evidence for a value 7 m i n ~ 50 in 
the jet of PKS 0637-752 based on spectral fitting the jet 
X-ray spectrum. We therefore adopt 7 m j n = 50 in the 
following analysis. 



The model and assumptions: — We use the standard one- 
zone IC/CMB model equations fsee iWorrall 120091 ) as- 
suming a continuous jet geometry (i.e. S v oc S 2+a for 
the synchrotron flux density). Many of the jet knots ap- 
pear elongated along the jet axis. Therefore, since this 
model requires the jet viewing angle to be small, we take 
account of projection effects in the calculation of knot 
volumes. Without independent constraints on the jet 
viewing angle, we simply assume a representative value 
in order to make an approximate de-projection. As we 
show in the following sections, the IC/CMB model re- 
quires a Doppler factor of order 6 > 6 at least in some 
parts of the jet. This implies that the jet viewing an- 
gle must be < 9 degrees. Angles significantly less than 
9 degrees are unfavourable since that would imply un- 



comfortably large source size. For example, if the jet lies 
closer than 5 degrees to the line of sight, the total source 
size must be > 2.3 Mpc. Therefore, we assume a jet 
viewing angle of 9 = 9 degrees in making an approximate 
de-projection of the knot length along the jet axis, and 
note that this angle corresponds approximately to the 
maximum jet viewing angle given the derived Doppler 
factor S > 6. The de-projected volumes, V, are related 
to the projected volumes Vb as V = VqJ sin#. The pro- 
jected volumes Vb are listed in Table [21 For the purposes 
of the current calculations, and to reduce the number 
of model parameters, we further assume that the proton 
contribution to the particle energy density is negligible. 

The results of IC / CMB modeling are presented in Ta- 
ble |U Included in the table is the jet kinetic energy 
flux, which we calculate based on the derived jet param- 
eters in the case of purely leptonic and electron/ proton 
j ets, us ing equation B17 of Appendix B in lSchwartz et al.l 
(|2006a|) . with the following simplifying assumptions: The 
bulk Lorentz factor of the jet T >> 1, equipartition be- 
tween particle and magnetic field energy densities, and 
a tangled magnetic field geometry such that (B' 2 ) = 
(2/3)B 2 . With these assumptions, the expression for 
kinetic luminosity in the case of a purely leptonic and 
electron/proton jet respectively, are (in c.g.s. units) 



LfJ~^7rR 2 T 2 c 



B* 
3- 

r- i 



ergs/s 



. n e m p c + — ergs/s 

1 / 07T ' 



In the case of an electron/proton jet, we assume one cold 
proton for each relativistic electron. The electron density 
n e , in c.g.s. units, is calculated from the spectral fits as 



B 2 



87rm e c 2 \ a — 1 



Train 



l 



where the magnetic field strength B is in Gauss. 

The model-dependent results presented in Table 0] sug- 
gest that, if the results of this model are correct, while 
the jet Lorentz factor remains approximately constant 
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along the jet, the magnetic field and particle density de- 
crease by a factor of a few between the innermost knots 
and the outer knots. The results indicate that there is 
modest, if any, loss of kinetic luminosity along the jet. 

3.3. Comparison of jet morphology with PKS 0637-752 
and PKS 0920-397 

The jet of PKS 2101-490 undergoes an apparent bend 
of approximately 45° at about half way between the 
core and jet termination, somew hat reminiscent of th e 
jet morphology in PKS 0637-752 (jSchwartz et al.ll2000f ). 
However, unlike PKS 0637-752, the X-ray emission in 
PKS 2101-490 is detected after the jet bend. This is of 
interest because a change in jet viewing angle produces 
a change in Doppler factor, which should manifest itself 
as change in knot brightness and the X-ray to radio flux 
density ratio. However, there is not a unique relationship 
between apparent bend angle and change in jet viewing 
angle. Therefore, the fact that X-ray emission continues 
beyond the bend in 2101-490 cannot be used to constrain 
the emission mechanism. It may simply be the case that 
the bend in 0637-752 is associated with a significant in- 
crease in viewing angle (and a hence significant decrease 
in the Doppler factor) whilst the bend in 2101-490 may 
be associated with a relatively small change, or a decrease 
in viewing angle (and hence small change, or increase in 
Doppler factor). Having said that, on a population basis, 
the probability distribution of change in jet viewing an- 
gle for a given apparent bend angle may enable a useful 
constraint on the emission mechanism. 

In addition to the morphological similarities to 
PKS 0637-752, we note a striki ng similarity betwee n 
PKS 2101-490 and PKS 0920-397 (ISchwartz et al.ll2010h . 
particularly in the vicinity of the jet termination. In both 
these sources, a bright, compact knot approximately 1 
arcsecond (> 10 kpc de-projected) upstream from the 
jet termination is detected in radio, optical and X-ray 
bands. A similar, X-ray bright pre-hots pot jet knot is 
seen m 1354+195 (jSambruna et al.ll2002D . One possible 
interpretation is that these pre-hotspot jet knots may 
be associated with an oblique shock that is produced as 
the jet enters the highly turbulent region at the head 
of the cocoon. Numerical simulations of extragalactic 
jets show that as the jet approaches the hotspot, it en- 
counters an increasingly violent environment, with strong 
turbulence that can perturb the jet flow, inducing insta- 
bilities or directly causing oblique shocks to f orm due 
to de nsity/pressure gradients in the lobes (eg. iNormanl 
1996). The pre-hotspot jet knots may be the result of 
the jet entering an increasingly violent environment near 
the hotspot. 

4. COMPARING INDEPENDENT ESTIMATES OF JET 
ENERGY FLUX FOR PKS 2101-490 

In this section, as a means to assess the validity of the 
IC/CMB model, our goal is to obtain independent esti- 
mates of jet power and compare these with the estimate 
of jet power derived from the IC/CMB model. 

4.1. Jet Energy Flux from Hotspot Parameters 

It is possible to estimate the jet power based on ob- 
served hotspot properties by applying the conservation of 
momentum and making a number of reasonable assump- 
tions about the hotspot. In the following section we carry 



out this approach, and compare the hotspot derived jet 
power estimates to those obtained from IC/CMB mod- 
elling of Knot 6 — the jet feature closest to the hotspot. 

Consider a uniform jet of area A, particle energy den- 
sity e, pressure p, mass density p, relativistic enthalpy 
w = e + p + pc 2 , magnetic field perpendicular to the flow 
direction B± (as indicated by the polarization), speed 
c/3 and corresponding bulk Lorentz factor T. The kinetic 
power (Ljct) and m omentum flux (Fm) along the jet are, 
in c.g.s. units (e.g. lDouble et al.ll2004T i: 



L jct = AT 2 l3c[w + 



B 2 



F, 



M l 



A 



T 2 /3 2 [ w + 



47T 

B 2 



4tt 



P 



Bl 
8tt 



(2) 
(3) 



Let us consider first the region of the jet upstream of 
the hotspot. In this region the IC/CMB model indicates 
that the jet Lorentz factor, T » 1. In such a highly 
relativistic jet the kinetic power £*E,jct is simply related 
to the momentum flux i^Mjet via 



L; ot sa c x Fm 



(4) 



We assume a near normal shock at the jet terminus and 
appeal to conservation of momentum, so that Fyijct — 
-FVhotspot, and thus L jct « c F Mi hots P ot- Assuming that 
the jet plasma in the hotspot is decelerated to a low 
velocity so that T 2 /3 2 <C 1, the following relation then 
holds: 



■cA x 



p 



Bl 
8tt 



(5) 



hs 



The above equality holds regardless of assumptions about 
the jet characteristics such as its composition or the ra- 
tio of magnetic to particle energy densities in the jet or 
hotspot. We estimate the jet kinetic luminosity simply 
by estimating the hotspot parameters. This technique of 
jet energy flux estimation will be considered further in 
an upcoming paper (Godfrey & Shabala, in prep.). 

In order to estimate hot spot parameters, we assume 
that the lepton population is ultra-relativistic (p = e/3) 
and dominates the particle pressure. Then 



Ac- 



B 2 



1 



1 e e ± 
3 e B 



(6) 



As a check on this analysis, we apply it to the case 
of Cygnus A. Synchrotron self Compton modelling of 
Cygnus A hotspot A, assuming a power law electron 
energy distribution of the form AT (7) — K e ~f~ a , indi- 
cates B=150 p, G and B Rn = 280 uG assuming R=2kpc 
and a = 2.1 (|Wilson et all 120001 ). For hotspot D: 
B=150 uG and B Rn = 250 /xG assuming R=2.2kpc and 
a = 2.05. iWilson et al.l (|2000f ) do not not give an es- 
timate of the electron energy density, so we estimate 
the electron energy density from the SSC and equiparti- 

b 2 ( r \( a+1 )/2 

tion magnetic field strengths as e e = \b~) » 
since e e oc K e and the synchrotron flux density S v oc 
K e B < - a + 1 ^ 2 =>• e e cx S-(° +1 )/ 2 , in the case of a power 
law distribution of the form AT (7) = K e ^ a . Therefore 
e e = 8 x 10~ 9 ergs/cm 3 and e e = 5.4 x 10~ 9 ergs/cm 3 
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TABLE 5 

Jet Power derived from IC/CMB modelling of 
K6 SED 



Composition 


Tmin - 


= 10 


7min = 


= 50 


Leptonic 


1 x 10 46 


ergs/s 


7 X 10 45 


ergs / s 


Electron /Proton 


3 x 10 47 


ergs/s 


4 X 10 46 


ergs / s 



for hotspots A and D respectively. Hence, we find 
Ljct ~lx 10 46 ergs/s for both hotspots A and D. This es- 
timate of the jet power in Cygnus A is in excellen t agree- 
ment with the value obtained by iWilson et al.l ([2006) , 
Ljct > 1.2 x 10 46 ergs/s, based on modeling the co- 
coon dynamics. An ind ependent estima te of jet power 
in Cygnus A comes from lLobanovl ([19981 1 who show that 
frequency dependent shifts of the radio core enable a de- 
termination of the jet parameters, which, when applied 
to the case of Cygnus A, gives a value for the jet power 
L jct = (0.55 ± 0.05) x 10 4b ergs/s. 

We now return to the case of PKS 2101-490. The 
Chandra X-ray image provides only an upper limit to the 
hotspot X-ray flux density. Therefore, to calculate the 
hotspot magnetic field strength and particle energy den- 
sity, we adopt equipartition estimates. We note that the 
adoption of an equipartition magnetic field strength does 
not overproduce hotspot X-ray emission, for which we 
have obtained an upper limit. We justify the assumption 
of equipartition by the fact that high luminosity hotspots 
typically exhibit X-ray flux density consistent with the 
equip artition synchrotron se lf Compton model predic- 
tions (jHardcastle et al.l 12004(1 . Moreover, the total en- 
ergy density, and hence derived jet power, is only weakly 
dependent on magnetic field strength for magnetic fields 
near equipartition conditions, so that moderate depar- 
ture from equipartition does not alter our conclusions. 
For example, if the hotspot magnetic field strength were 
1/3 the equipartition value, then our equipartition jet 
power would under-estimate the true value by less than 
a factor of 2. 

The hotspot synchrotron spectral index is not well con- 
strained, so we apply a single power law model along 
with the following assumptions: 10 < 7 nml < 1000; 
7max > 10 4 ; a = 0.8. This model, combined with 
the flux densities and volume listed in Table [2] indi- 
cates that -B eq ,hs ~ 200 — 450 (j,G. Due to the prior 
assumption of equipartition we have e e ± = ee- Taking 
A = (6 ± 1) X 10 43 cm 2 as estimated from the 20 GHz 
ATCA image, the jet energy flux calculated from the 
hotspot model is 



L jct = (0.4 - 2) x 10 4b ergs/s 



(7) 



The range of values reflects the range of 7 nml assumed in 
the hotspot model, with the lower jet power correspond- 
ing to higher assumed values of 7 m i n - Minimum Lorentz 
factors 7r nin ~ 600 have been observed in a number of 
hotspots ([Godfrey et al.1 [20091 and references therein). 
This is an order of magnitude greater th an the value of 
7m in a ssumed in the jet. As suggested by lGodfrev et al.l 
(2009), such an increase in 7 m i n between the jet and 
hotspot may be due to the dissipation of jet energy in 
a relativistic proton/electron jet terminating at a near- 
normal shock. 



4.2. Willott et al. Q ]c t — £151 relation 

A widely used method for estimating jet power in 
high luminosity rad io sources is based on the model of 
IWillott et ail (| 1999ft . in which the jet power is estimated 
using the 151 MHz radio luminosity as follows: 



Qw « / 3/2 1.5 x 10 38 ( 



10 28 W HZ- 1 sr~ 



6/7 



W (8) 



where Qw is the kinetic power per jet, Si 51 is the flux 
density of the entire source at 151 MHz, Dl is the lu- 
minosity distance, and f is a parameter that accounts 
for uncertainties in various model assumptions, with 
1 < / £ 20 (see IWillott et all Il99l for details of the 
model). For high power sources, the value of f has not 
been constrained empirically, but is typically taken to be 
/ = 1 - 20 (e.g. iFernandes et alj|2011t jHardcastle et al.l 
12007(1 . We extrapolat e flux density me asurements of 
S408 MHz = 2.34 Jy (lLarge et al.l 119811 ). S 8 43 mh z = 
1.48 Jv (iMauch et all 12003ft and S1.41 G h z = 1.1 Jy 
([Wright fc Otrupceld 119901 ) to estimate the flux density 
at 151 MHz — S151 mhz ~ 4.25 Jy — and thereby obtain 
Qw - 1 (3) x 10 46 ergs/s assuming / = 10 (20). 

4.3. Average jet power from lobe parameters and source 

age 

Here we estimate the average jet power over the life of 
the source based on the lobe parameters and an estimate 
of the source age. 



i-agc 

where U\ h e is the total (magnetic plus particle) energy 
contained in the lobe. The factor of 2 is an approx- 
imate correction accounting for the fraction of jet en- 
ergy converted to lob e kinetic energy and work done by 
the e xpanding lobes ([Willott et al.l 119991 : iBicknell et al.l 
11997( 1. In SjTjwe estimate the lobe parameters via syn- 
chrotron and inverse Compton modelling of the radio 
to X-ray SED, and argue that the source age is likely 
^agc ^ 10 yrs. Based on the results of $7] we estimate 
the average jet power over the life of the source to be 
-kjet si 10 45 ergs/s. This is a factor of 4 - 20 lower than 
the other estimates of jet power. We note that we have 
ignored the possible contribution of protons to the en- 
ergy density of the lobes, and we have ignored adiabatic 
cooling in estimating the cooling rate in the above analy- 
sis, each of which would increase the jet power estimate. 
We note that the calculations in $7] involve projected vol- 
umes. Using de-projected volumes will not increase the 
total total energy, since the product K e x V remains 
constant, and the particle energy dominates Ui o h in this 
source. 

4.4. Discussion of the Results 

Three independent order-of-magnitude estimates of jet 
power indicate Lj Ct ~ 10 45 — 10 46 ergs/s. For comparison, 
the jet power obtained from the one-zone IC/CMB model 
of jet X-ray emission are given in Table [5j for different 
values of the minimum Lorentz factor "f m in, and different 
assumed jet compositions. 

It can be seen that, in the case of a purely leptonic com- 
position, there is good agreement between the IC/CMB 
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Fig. 7. — Spectral energy distributions for knots 3 (top), 4 (middle) and 6 (bottom). Also plotted are the synchrotron and analytic power 
law model for the IC/CMB interpretation, with a = 0.8, 7 m i n = 50, and 7 ma x = 10 , and in the case of knot 6, a broken power law model 
is used (see j5]for details). The parameters for the fits are given in Table [4] Note that while the X-ray spectral index of individual knots 
could not be measured, the X-ray spectral index of the entire jet is j e ^ v = 0.85 ± 0.2, consistent with the slope of the model in the 
X-ray band. 

10 46 ergs/s for both leptonic and electron/cold proton 
compositions, provided 7 m i„ > 50. An important point 
to note about our jet power estimates is that, for the 
IC / CMB model, we have used a cylindrical geometry for 
knot K6 with length approximately 10 times its width. 
This highly elongated cylindrical geometry is required 
by the observation that knot K6 is elongated along the 
jet axis in the radio maps, and the fact that the jet 
must be aligned within < 9 degrees of the line of sight 
if the IC/CMB model is to be applied. Typically, jet 
knots are assumed to be spherica l with the knot diam- 
eters equal to the jet widths (e.g. iTavecchio et"aT1l2000l : 
IMehta et all 12009 s ). If we were to assume a spherical 
volume in our analysis, the IC/CMB derived jet powers 
would be greater by a factor of ~ 4. 

5. THE OPTICAL EMISSION FROM KNOT 6 

5.1. Introduction 

The rate of optical jet discoverjO in objects with con- 
firmed radio and X-ray jets is approximately 2 out of 3. 
However, few of these sources have resolved optical knots, 
particularly at high redshift (excepting P KS 0637-752 in 
which multiple optical kn ots are resolved: iSchwartz et al.l 
120001 IMehta et al.ll2009h . and whilst a systematic study 
has not yet been performed, it would appear that rela- 
tively few have detailed information regarding the spec- 
tral slopes in both the radio and optical bands. 

In PKS 2101-490, knot 6 is the only jet region that is 
reliably detected in all three bands (radio, optical and X- 
ray). The radio through optical spectrum for knot 6 can 
be fit using a broken power law with a break of Aa = 0.5 



10" 10" 10" 10" 

Frequency [Hz] 



Fig. 8. — Broken power law synchrotron spectrum fit to the radio 
to optical spectral energy distribution of knot 6. The relevant 
parameters of the fit are: a = 0.8, v\y = 7 X 10 12 Hz, i/ max = 
10 16 Hz. 

model and independent estimates of jet power, regardless 
of the value for 7 m i n within the stated limits. Under the 
assumption of a cold proton/electron composition, there 
is good agreement between the IC/CMB model and in- 
dependent estimates of jet power, provided jxain ^ 50. 
This assumes one cold proton per radiating electron. If 
relativistic protons are included, the lower limit on 7 m i n 
would increase. 

We have shown that in the case of PKS 2101-490, the 
IC/CMB predicted jet power is << 10 48 ergs/s (the "un- 
comfortably large" jet power derived for some quasar X- 
ray jets based on the IC/ CMB model (|Atovan fc Dermerl 
120041: IMehta et all 12009( 1) under a range of reasonable 
assumptions. Indeed, the IC/CMB derived jet power is 
comparable to independent estimates of Lj ~ 10 45 — 



17 See http:/ /astro. fit.edu/jets/ for a listing of known optical/IR, 
jet sources, and information tor each source. 
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(Figure [7]). This suggests that radiative cooling may be 
responsible for the steeper spectrum at optical frequen- 
cies. However, the interpretation of the optical emission 
from Knot 6 is complicated by fact that the optical to X- 
ray spectral index is 32xio 14 Hz ~ ^-8, consistent with 
the observed radio and X-ray spectral index of the entire 
jet. It is therefore plausible that a significant fraction of 
the observed optical flux density is IC/CMB emission. 
For the IC/CMB model to produce the observed optical 
flux at 6.32 x 10 14 Hz, the electron energy distribution 
must continue with the same slope to Lorentz factors 
7 < 10, assuming r ~ S = 6. If the electron distribu- 
tion cuts off at a Lorentz factor 7 >> 10, the IC/CMB 
emission will not make a significant contribution to the 
observed optical flux density. 

With this caveat in mind, in this section we model 
the radio through optical spectrum in terms of a bro- 
ken power law synchrotron model, and consider whether 
the jet parameters derived from IC/CMB modeling are 
consistent with this interpretation of the radio to optical 
spectrum of knot 6. Figure [3] illustrates the broken power 
law fit to the SED of knot 6. The parameters of this fit 
are a = 0.8, v h = 7 x 10 12 Hz and f max = 10 16 Hz. 

5.2. Synchrotron model for the radio to optical spectrum 

of Knot 6 

In modeling the radio to optical s pectrum we consider 
the e ffects of synchrotron cooling (Mciscnhcime r et al.l 
1989) and adopt the following procedure: We assume 
that relativistic electrons are injected at a shock with an 
energy index a (number density per unit Lorentz factor, 
AT (7) = A" e 7~ a ). The electrons cool as a result of syn- 
chrotron and inverse Compton emission downstream of 
the shock over a length D, which defines the extent of the 
emitting region. In this model we assume that the post- 
shock magnetic field and velocity are uniform; this is a 
reasonable assumption for an oblique shock. An electron 
with Lorentz factor 7 cools as a result of synchrotron and 
inverse Compton emission according to: 



where 



dt 

where t' is time in the plasma co-moving frame, 

4 u T ( B 2 TT \ 
\ — + U cmb J , 



3 m e c \ 8tt 



(10) 



(11) 



in c.g.s. units, and the energy density of the cosmic mi- 
crowave background in the plasma co-moving frame is 
given by 

Ucmb = Ucmbv(1 + z) a T 2 (12) 



CMB energy density is 
~ 3 and r is the jet Lorentz 



where the current epoch 
UcuBfl ~ 4.2 x 10~ 13 ergs cm 
factor. 

The integrated emission from the post-shock region is 
the superposition of the emission from a number of pro- 
gressively cooled slices. The volume averaged number 
density per unit Lorentz factor is described by: 





N( 7 ) = I K Blb . 

(a-l) 



-(«+!), 



(i) 



7 < 7min, 7 > In 
7min < 7 < 7max 



7 < 76 
7 > 7ft 



(14) 



The break Lorentz factor 7^ is given by: 

7ft = ^r sh /3 sh (15) 

where /3 s h is the velocity of plasma relative to the shock 
and T s h is the corresponding Lorentz factor. The electron 
Lorentz factor 7b is the Lorentz factor to which an elec- 
tron of initially infinite energy cools following the shock. 
If the jet-shock is stationary then T s h = T. 

The electron distribution (I13[) describes a broken power 
law spectrum with the electron spectral index smoothly 
changing from a to (a + 1) at 7 w 7^. The correspond- 
ing synchrotron spectrum is a broken power law with 
spectral index smoothly changing from a to a + 0.5 at a 
frequency in the observer's frame, 



(1 + Z) 47T 



(16) 



where r2o(= m c.g.s. units) is the non-relativistic 
gyro-frequency. This equation relates the break fre- 
quency to the break Lorentz factor, and it is the break 
frequency, i>b, that is actually used in the fit; the other 
parameter is the injected spectral index, a — (a — l)/2. 
The values for the magnetic field B ss 30 /iG and Dopplcr 
factor 6 ss 6 have previously been determined from the 
relationship between X-ray and radio emission in § 13.2.21 
In order to minimize the number of parameters we as- 
sume that r = (5 in the following analysis. 

We can compare the fitted break frequency with the 
value theoretically implied by equations (|15|) for 7^ and 
(|16[) for V]j. The theoretical break frequency is 



27 



Vb = 



64tt 1 + z D 2 a 2 



F 2 B 2 

1 sh^sh 
U 2 



(17) 



(13) 



where U = U-q + C^cmb is the total magnetic plus radi- 
ation energy density. The jet parameters derived from 
the IC/CMB model for this knot imply that inverse 
Compton losses dominate over synchrotron losses, since 
the CMB energy density in the jet rest frame is greater 
than the magnetic energy density. In order to compare 
the predicted break frequency with that determined by 
the synchrotron model we need to estimate the length, 
D, of the emitting region. The radio image of knot 6 
is found to be elongated in the jet direction with de- 
convolved (intrinsic) length D sa 0'.'6. The upper limit 
to the size of the X-ray knot is 0'.'6. The maximum 
allowed viewing angle for 8 = 6 is 9 max s» 9°, there- 
fore, we take the de-projected knot length to be of order 
D w 0"6/sin6> max w 4" « 30 kpc « 9 x 10 22 cm. If 
we adopt D s» 10 23 cm and insert in equation (fT?]) the 
appropriate parameter values from the IC/CMB fit we 
obtain a predicted break frequency of 1.2 x 10 13 Hz. This 
is less than a factor of two above the modelled value 
of 7 x 10 12 Hz. Hence, the jet parameters derived from 
IC/CMB modeling of knot 6 are broadly consistent with 
the break frequency estimated by fitting a broken power 
law to the radio and optical data points. 
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We note that the F814W HST image of knot 6 (Figure 
[3]) contains a double-peaked appearance, suggesting that 
a one-zone model of the emission region, which spans 
>30 kpc in the radio band, may not be applicable. 




2 3 4 5 

Energy (keV) 



Fig. 9.— Observed ACIS-S spectrum of the core of PKS 2101-490 
with model spectrum (absorbed power-law) overlaid. Also shown 
are the residuals of the fit, expressed as [(data - model) /error] in 
each bin. The positive residuals at approximately 3.1 keV may be 
associated with a 6.4 keV Fe Ka lin e redshifted by z=1.0 4, consis- 
tent with the spectroscopic redshift (Marshall et al. 2005, Gelbord 
& Marshall, in prep.). 

6. X-RAY SPECTRAL ANALYSIS OF THE QUASAR CORE 

X-ray events were extracted using a circular aperture 
of radius 2" centred on the core peak position. A to- 
tal of 3687 events were extracted with energies in the 
range 0.5-7 keV. The pile-up fraction is estimated using 
the PIMMS Proposal Planning Toolkit to be 2%, and is 
therefore neglected in the following analysis. The pulse- 
height amplitude spectrum was re-binned to have a min- 
imum of 20 counts per bin, so that the % 2 statistic could 
be applied. 

The spectrum was fit using an absorbed power law 
model, with HI column density fixed at the Galactic 
value A^hi = 3.41 x 10 20 cm -2 as determined from the 
COLDEN column density calculator provided by the 
Chandra X-ray Center. The data, along with the best 
fit curve are shown in Figure [9] The parameters of the 
best fit are a X -ra y = 0.84 ± 0.05, S = 87 ± 3 nJy at 
1 keV, and chi-squared per degree of freedom \ 2 / d.o.f. 
= 146.4 / 147 ^ 1 00. 

iMarshall et~al"1 (|2005D reported a redshift of z = 1.04 
for this source, based on an unpublished Magellan spec- 
trum (see JQ) ■ We note positive residuals consistent with 
a 6.4 keV Fe Ka line redshifted by z ss 1.04. If we in- 
clude a line in the chi-squared fit, the model parameters 
become a X -ra y = 0.85 ± 0.05, S = 87 ± 3 nJy at 1 keV, 
and the chi-squared per degree of freedom becomes % 2 / 
d.o.f. = 134.7 / 145 « 0.93. Alth ough an F-test would in - 
dicate that the line is significant. iProtassov et all ()2002f ) 
show that such a test is not applicable to adding a spec- 
tral line to our basic model, and we have not carried out 
the extensive Monte Carlo testing they recommend since 
the jet physics is the primary result of this paper. 

7. LOBES 



Within the lobes, the radiation energy density of the 
CMB (Uqmb — 7 x 10 -12 ergs cm -3 ) is greater than the 
energy density of the locally generated synchrotron radi- 
ation. We therefore model the X-ray emission from the 
radio lobes in terms of inverse Compton scattering of the 
CMB, in this case assuming 6=1. In modeling the radio 
lobe emission we assume a broken power law electron en- 
ergy distribution of the form given by equation 1131 The 
value for the break Lorentz factor is constrained by the 
X-ray spectrum. We argue as follows. The radio spectral 
index of the lobe is a r = 1.05 ±0.1. We therefore assume 
that the radio emission corresponds to the post-break re- 
gion of the electron spectrum. We assume a low-energy 
electron index, a = 2.1, so that the slope of the model 
synchrotron spectrum above the break matches the ob- 
served radio spectral index. The X-ray spectral index in 
the lobe is ax-ray = 1.3 ± 0.3. The inverse Compton 
X-ray emission in the Chandra X-ray band (0.5 - 7 keV) 
is produced by 7 « 600 — 2000 electrons, given that the 
scattered photons are CMB photons with an average en- 
ergy fa 3.6 7 2 kT C MB (1 + z), where Tomb ~ 2.7 K is 
the current epoch temperature of the cosmic microwave 
background. The steep X-ray spectrum therefore indi- 
cates that the Lorentz factor corresponding to the cool- 
ing break must be 7^ < 600. Hence, we adopt 7b = 500. 
Given the unknown jet viewing angle, we assume pro- 
jected volumes in the analysis. If the jet is angled close 
to the line of sight, as required by the IC/CMB model for 
jet X-ray emission, the volumes may be underestimated 
by a factor of 5 - 10. Larger volumes would imply closer 
agreement with the equipartition model. 

Figure [TU] shows the model fit to the radio and X-ray 
data for the lobe and counter lobe. The best-fit parame- 
ter values are given in Table [51 along with the minimum 
energy magnetic field strength calculated usi ng standard 
expressions (e.g. Worrall & Birkinshaw 2006). The mag- 
nitude of the average magnetic field strength, as well as 
the ratio of B to B m j n is typical of powerful radio galaxies 
(jCroston et al.ll2005f) . 

The presence of a low value for the break Lorentz factor 
(7b ^ 500) indicates an old source. However, enhanced 
inverse Compton cooling leads to a shorter cooling time 
than for radio galaxies at the current epoch. The cooling 
time for electrons of Lorentz factor 7 in a magnetic field 
B (and hence magnetic energy density Ub) and immersed 
in the microwave background radiation with energy den- 
sity U C MB is 

m e c x f u B y 1 , 
^• 6xl ° 8 ( 1 + Z^)"(5oo)"^^ 

for a microwave background temperature (1+z) x 2.725K. 
For the estimated lobe magnetic field strength, Ub ~ 
C^cmb, and 76 ~ 500, implying a source age t ag0 < 
2 x 10 8 yrs. Adiabatic losses would increase the cool- 
ing rate, meaning the above prediction may be an over- 
estimate of the source age, hence the inequality in the 
estimate of £ age given above. For a deprojected jet length 
> 500 kpc (viewing angle < 9° as required by the 
beamed, equipartition IC/CMB model for jet X-ray emis- 
sion) this implies a source expansion speed > 0.01 c. This 
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Fig. 10. — Lobe and counter-lobe spectral energy distributions. The solid lines are the best fit synchrotron plus IC/CMB model. The 
parameters of the best fit models are given in table \E\ The "bow-tie" around the counter-lobe X-ray data point illustrates the constraint 
on the X-ray spectral index of a x = 1.3 ± 0.3. There are not enough X-ray counts in the lobe to provide meaningful constraints on the 
spectral index. However, in modelling the lobe SED we have assumed the same X-ray spectral index as the counter-lobe. 



TABLE 6 

Parameters of model fits to the lobe spectra. 



Volume 


a 


Tin in Tb 


Tin ax 


K e 


B Smin 


[cm 3 ] 








[cm -3 ] 





Counter-lobe 9 x 10 69 0.55 10 500 6 x 10 5 1 x 10~ 5 14 29 
Lobe 4 x 10 69 0.55 10 500 6 x 10 5 0.8 x 10~ 5 13 25 



estimate does not differ greatly from the estimate of the 
typical expans ion speeds of l obes of a few percent of the 
speed of light (|Scheuerlll995[ ). 

8. CONCLUSIONS 

We have presented an analysis of new Chandra, HST, 
and ATCA images for the quasar jet source PKS 2101- 
490. We extracted the radio to X-ray spectral energy dis- 
tributions from seven regions of the 13" jet, and modeled 
the jet X-ray emission in terms of Doppler beamed in- 
verse Compton scattering of the cosmic microwave back- 
ground (IC/CMB) for a jet in a state of equipartition be- 
tween magnetic and particle energy densities (33|. Using 
this model, we derived a bulk Lorentz factor T ~ 6 and 
magnetic field strength of order 30 /iG. 

A major goal of this work has been to assess the valid- 
ity of the beamed, equipartition IC/CMB model for jet 
X-ray emission through the use of additional constraints: 
namely, independent estimates of jet power, and the lo- 
cation of an apparent cooling break in the synchrotron 
spectrum of one of the jet knots. Below we discuss the 
results of our analysis for each of these additional con- 
straints. 

The jet power predicted by the Doppler beamed, 
equipartition IC/CMB model was found to be in good 
agreement with independent order-of-magnitude esti- 
mates of jet power for this object, provided that 7 m ; n ~ 
50 in the jet, and the knots are significantly longer than 
the jet width, as implied by de-projection of the observed 
knot lengths (gj. 

The brightest X-ray knot was detected in two HST 
filters, and the radio to optical data points were modelled 



as a broken power law with a standard cooling break 
Aa = 0.5. The inferred break frequency was found to be 
consistent with the break frequency predicted using the 
IC/CMB model parameters along with the de- projected 
knot length that is implied by the small jet viewing angle 
required by the IC/CMB model O- However, we noted 
sub-structure in the F814W HST image of the bright 
optical knot, which is incompatible with a simple one- 
zone continuous injection model. 

Finally, we note positive residuals consistent with a 
6.4 keV Fe Kg line redshif t ed by the estimated z=1.04 
reported bv lMarshall et all (|2005[ ). 
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